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ABSTRACT 


Vhis  report  presents  an  interactive  Fortran  77  computer  program  for 
developing  and  plotting  standard  Shewhart  control  charts  for  X-bar  and  R 
or  s,  p,  np,  c  or  u.  It  was  developed  on  a  VAX  11/750  computer  using  the 
VMS  operating  system.  Central  lines  and  control  limits  may  be  based  on 
the  input  data,  modifications  of  input  data,  or  standard  values.  Input 
data  is  stored  in  separate  named  files  so  that  it  may  be  recalled,  modified 
or  added  to.  * 


USERS  GUIDE 


Before  Running  Program 

All  data  must  be  stored  in  a  permanent  data  file.  The  name  of  the  data 
file  is  limited  to  10  characters. 

X-bar,  R,  and  S  charts 

This  data  must  be  stored  in  a  data  file  with  one  subgroup  per  line  with  a 
space  or  a  comma  between  each  value.  The  maximum  number  of  subgroups  is  30, 
and  the  maximum  subgroup  size  is  10.  The  maximum  data  value  permitted  is 
9999.99. 

P  Charts 


This  data  must  be  stored  in  a  data  file  in  the  order  number  inspected, 
number  rejected,  separated  by  a  comma  or  a  space.  The  maximum  number  of  data 
pairs  is  30.  The  maximum  data  value  permitted  is  9999.99. 

U  Charts 


This  data  must  be  stored  in  a  data  file  in  the  order  number  of 
nonconformities,  number  produced,  separated  by  a  comma  or  a  space.  The 
maximum  data  value  permitted  is  9999.99. 

NP  and  C  Charts 


This  data  must  be  stored  in  a  data  file  with  data  values  separated  by  a 
space  or  a  comma.  The  maximum  number  of  data  points  is  30.  The  maximum  data 
value  permitted  is  9999.99. 

If  previously  used  control  limits  are  going  to  be  used,  then  they  must  be 
stored  in  a  permanent  data  file.  They  must  be  stored  in  the  order-upper 
control  limit,  lower  control  limit,  central  line,  separated  by  a  space  or  a 
comma.  If  both  X-bar  and  R  (S)  control  limits  are  to  be  used,  they  must  be 
stored  in  separate  files. 

Running  the  Program 

There  are  2  basic  menus  in  this  program.  The  first  menu  is  only  shown 
once.  It  gives  the  user  a  choice  of  charts.  If  the  user  should  wish  to 
change  the  type  of  chart,  the  program  will  have  to  be  run  again.  This  menu 
gives  the  user  a  choice  between:  1)  X-bar  and  R  charts,  2)  X-bar  and  S  charts, 
3)  C  or  U  charts,  4)  P  or  NP  charts.  The  user  makes  his/her  choice  by 
entering  the  number  associated  with  the  option  and  hitting  the  carriage 
return. 

After  this  choice  has  been  made,  the  user  will  have  to  enter  the  names  of 
the  data  file  and  output  file.  These  names  are  limited  to  10  characters. 


X-bar,  R,  S  Charts 


The  user  must  enter  the  number  of  subgroups  and  the  subgroup  size.  These 
values  must  be  limited  to  30  and  10,  respectively.  The  name  of  the  X-bar  and 
R  (S)  control  limit  files  must  also  be  supplied.  Even  if  these  files  do  not 
exist  at  this  time,  the  user  is  still  required  to  enter  a  name  for  them  in 
case,  during  the  run  of  this  program,  a  decision  is  made  to  save  current 
control  limits. 

P  and  NP  Charts 


If  the  user  chooses  this  as  the  first  option,  he/she  must  then  enter  a  1 
if  working  with  fraction  rejected  or  2  if  working  with  number  rejected.  The 
user  must  also  give  the  name  of  the  control  limit  file  (limit  10  characaters) . 
The  number  of  lots  inspected  must  be  entered  and,  if  using  the  NP  chart,  the 
lot  size  must  also  be  given. 

C  and  U  Charts 


If  the  user  chooses  this  as  the  first  option,  he/she  must  then  enter  a  1 
if  working  with  number  of  nonconformities,  or  a  2  if  working  with 
nonconformities  per  unit.  The  number  of  data  points  (pairs)  must  be  entered 
along  with  the  name  of  the  control  limit  file. 

All  of  the  above  information  (and  future  information)  will  be  requested 
by  the  computer  in  the  form  of  a  written  prompt  which  will  appear  on  the 
screen.  The  user  must  wait  for  the  prompt  before  entering  the  information. 

After  the  previous  information  has  been  entered,  a  second  menu  will 
appear  on  the  screen.  This  menu  will  vary  according  to  the  type  of  chart 
chosen  in  the  first  menu.  All  options  are  listed  below  plus  any  information 
that  the  user  will  need  to  enter. 

Options 

Hardcopy  of  Data 

A  copy  of  the  raw  data  and  the  points  plotted  on  the  charts  will  be 
written  to  the  output  file. 

Histogram 

This  option  will  only  appear  if  the  user  choose  X-bar  and  R  (S)  charts  in 
the  first  menu.  The  user  will  have  to  decide  if  the  histogram  will  appear  on 
the  screen,  written  to  the  output  file,  or  both.  The  number  of  cell  intervals 
will  also  have  to  be  entered.  If  the  histogram  is  viewed  on  the  screen,  the 
user  will  have  to  enter  a  1  to  continue  the  program. 

Control  Chart 


If  the  user  is  running  a  P  chart,  then  he/she  has  a  choice  between  using 
individual  lot  sizes,  estimated  average  lot  size,  or  calculated  average  lot 
size. 


iv 


1 


All  charts  give  the  op-ion  of  using  control  limits  calculated  from  data, 
old  control  limits,  or  aimed  at  values,  except  in  the  case  of  np  (individual 
lot  sizes)  charts  or  U  charts  which  do  not  give  the  option  for  using  old 
control  limits.  In  the  latter  case,  the  central  line  of  a  previous  control 
chart  could  be  entered  as  an  aimed  at  value.  If  aimed  at  values  are  chosen, 
these  must  be  entered  by  the  user. 

Charts  can  be  viewed  on  the  screen,  written  to  a  file  or  both.  The 
control  limits  can  also  be  saved  for  future  use.  Out  of  control  points  are 
indicated  by  a  question  mark. 

Add  More  Data 

The  user  must  enter  the  number  of  new  data  points  (or  subgroups).  This 
number  is  limited  to  30.  The  permanent  data  file  will  be  altered  to  reflect 
the  change  in  data.  If  the  total  number  of  data  points  (old  and  new)  is 
greater  than  30,  some  of  the  beginning  data  points  will  be  deleted  to  make 
room  for  the  new  ones. 

Delete  Data  Points 

This  option  allows  the  user  to  delete  data  points  from  the  calculation  of 
control  limits.  The  user  must  give  the  total  number  of  data  points  to  be 
deleted  and  the  number  of  each  point  to  be  deleted.  Once  this  is  done,  the 
control  chart  is  automatically  calculated,  so  the  user  will  have  the  same 
options  listed  under  Control  Chart.  These  deleted  points  will  be  plotted  on 
the  chart.  The  user  should  choose  to  have  the  control  charts  calculated  from 
data  or  else  the  influence  of  the  deleted  points  will  not  be  shown. 

Exit  Program 

This  option  ends  the  program.  To  obtain  a  hardcopy  of  everything  written 
to  the  output  file,  the  user  should  type  PRINT  (filename). 


INTRODUCTION 


The  following  program  is  an  interactive  program  designed  for  the 
development  of  X-bar  and  R  or  S,  P,  NP,  C,  and  U  charts.  The  program  features 
options  that  enable  the  user  to  analyze  data  using  standard  Shewhart  quality 
control  procedures.  Options  include  a  histogram,  the  ability  to  delete  points 
from  the  calculation  of  control  limits,  and  the  ability  to  add  new  points. 
Output  from  this  program  (excluding  hardcopy  of  the  data)  can  be  viewed  on  the 
screen  and  all  can  be  stored  in  a  data  file  to  be  printed  out  after  the 
program  run  is  finished.  This  program  is  written  in  Fortran  77  and  was 
implemented  on  a  Vax  11/750  computer.  It  is  intended  as  an  instructional  tool 
and  as  a  research  tool  for  analyzing  process  data. 


MAIN  PROGRAM 


DESCRIPTION  OF  PROGRAM 


The  basic  function  of  the  main  program  is  to  identify  the  type  of  data 
being  used,  read  the  data  from  an  external  data  file,  and  to  call  the 
appropriate  subroutines  depending  on  the  options  specified  by  the  user. 


The  first  menu  is  the  main  program  identifies  the  type  of  data  being 
used.  This  menu  will  only  be  displayed  once  during  a  run  of  the  program,  so 
if  the  user  wishes  to  change  the  type  of  data  being  used,  the  program  will 
have  to  be  run  again.  The  choices  offered  in  the  first  menu  are:  (1)  X-bar 
and  R  charts;  (2)  X-bar  and  S  charts;  (3)  C  or  U  charts;  (4)  P  or  NP  charts. 

Some  information  will  have  to  be  entered  by  the  user  from  the  terminal. 
This  information  includes  the  name  of  the  data  file,  the  name  of  the  output 
file,  the  name  of  the  control  limit  file  (except  in  the  case  of  a  u  chart), 
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the  subgroup  size  and  the  number  of  subgroups  (in  the  case  of  X-bar  and  R  or 
X-bar  and  S  charts),  and  the  number  of  data  points  (for  all  other  charts). 

The  data  will  then  be  read  from  the  external  data  file. 

The  second  menu  then  will  appear  on  the  screen  and  the  options  will  vary 
according  to  the  type  of  data  being  used.  Once  an  option  is  chosen, 
subroutines  corresponding  to  these  options  are  called  and  executed.  When  the 
option  is  completed,  this  second  menu  is  again  displayed. 

OPTIONS 

Hardcopy  of  data 

This  option  will  give  the  user  a  copy  of  the  current  data  set  and  the 
values  of  the  points  which  are  to  be  plotted  on  the  control  chart.  If  the 
data  type  being  used  corresponds  to  a  C  chart  or  an  NP  chart,  then  just  the 
subroutines  to  print  a  hardcopy  of  the  data  is  called.  If  any  other  type  of 
data  is  used,  then  the  first  subroutine  called  is  the  one  that  calculates  the 
values  of  the  points  that  will  be  plotted  on  the  control  chart.  This 
subroutine  then  returns  to  the  main  program  and  the  corresponding  subroutine 
to  print  a  hardcopy  of  the  data  is  called.  All  this  information  is  written  to 
the  output  file  named  by  the  user,  and  a  hardcopy  can  be  obtained  after  the 
user  has  finished  running  the  program. 

Histogram 

This  option  will  yield  a  histogram  either  printed  on  the  terminal  screen 
or  written  to  the  output  file  or  both.  This  option  is  only  available  if 
variables  data  are  being  used  to  prepare  X-bar  and  R  charts  or  X-bar  and  S 
charts.  When  this  option  is  completed,  the  second  menu  will  appear  on  the 
screen. 


2 


Control  Chert. 

In  general,  this  option  will  plot  a  control  chart  on  the  terminal  screen 
or  write  it  to  the  output  file,  or  both.  No  matter  which  control  chart  is 
requested,  three  subroutines  are  called  to  produce  the  control  chart.  First, 
a  subroutine  is  called  to  calculate  all  values  needed  to  compute  the  control 
limits.  These  values  are  then  passed  to  a  subroutine  that  calculates  the 
control  limits  and  these  values,  along  with  the  points  to  be  plotted,  are 
passed  to  the  subroutine  that  prints  out  the  charts.  When  this  option  is 
completed,  the  second  menu  again  will  appear  on  the  screen. 


Addina  new  subroutines/Adding  more  data 


These  options  enable  the  user  to  alter  the  data  currently  being  used.  A 
subroutine  corresponding  to  the  type  of  data  being  used  is  called.  The  user 
enters  the  new  data  from  the  terminal.  The  external  data  file  is  altered  to 
accomodate  the  new  data  as  is  the  current  data  array  in  the  program.  When 
this  option  is  completed,  the  second  menu  is  displayed  on  the  screen. 


Deleting  data  points 


This  option  enables  the  user  to  alter  the  values  of  the  control  limits  by 
eliminating  points  from  the  calculation  of  the  control  limits.  These  points 
are  NOT  deleted  from  the  current  data  array  in  the  program  or  from  the 
external  data  file.  The  first  subroutine  called  calculates  the  data  used  in 
the  calculation  of  the  control  limits  using  all  the  data  points.  The  next 
subroutine  called  alters  these  calculation  to  reflect  the  change  introduced  by 
the  deletion  of  points.  These  altered  values  are  then  passed  to  the 
subroutine  that  calculates  the  control  limits  and  the  points  to  be  plotted  are 
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then  passed  to  the  subroutine  that  plots  the  control  charts.  When  the  option 
is  finished,  the  second  menu  will  appear  on  the  screen. 


Exit  Program 

This  option  ends  the  program.  The  user  can  then  obtain  a  hardcopy  of  the 
output  file  which  will  contain  everything  that  was  written  to  the  file  (i.e. 
hardcopy  of  data,  histogram  and/or  control  charts  as  applicable  and  as 
requested). 

In  the  above  sections,  the  options  and  the  various  subroutines  connected 
with  them  were  described  in  general.  Each  subroutine  now  is  described  in 
detail. 


DESCRIPTION  OF  SUBROUTINES 
Subroutine  CRDXRS.  CROPU,  CRONPC 

These  three  subroutines  take  care  of  the  intermediate  calculations  needed 
to  compute  control  limits.  They  are  called  not  only  when  a  control  chart  is 
requested,  but  also  when  a  hardcopy  of  the  data  is  needed,  and  when  data 
points  are  to  be  deleted  from  control  limit  calculations. 

Subroutine  CROXRS 

This  subroutine  performs  calculations  on  subgroup  data.  The  raw  data, 
the  number  of  subgroups  and  the  subgroup  size  is  passed  to  the  subroutine. 

The  average  and  arange  (std.  dev.)  of  each  subgroup,  the  grand  average  (over 
all  subgroups),  the  average  range  (std. dev.)  are  passed  back  to  the  main 


program. 


Subroutine  CRDPU 


This  subroutine  is  called  to  perform  calculation  on  data  dealing  with  the 
number  rejected  and  number  inspected,  or  the  number  of  nonconformities  and  the 
number  produced.  The  raw  data  and  the  number  of  data  points  are  passed  to 
this  subroutine.  For  each  pair  of  data  points,  either  the  fraction  rejected 
or  the  number  of  nonconformities  per  unit  is  calculated.  The  average  fraction 
rejected  or  the  average  nonconformities  per  unit  is  calculated. 

If  the  data  being  used  deals  with  the  number  rejected  and  the  number 
inspected,  then  the  user  has  a  choice  of  using  individual  lot  sizes,  estimated 
average  lot  size,  or  a  calculated  average  lot  size.  If  the  user  chooses  a 
calculated  average  lot  size,  then  this  value  is  calculated.  If  the  user 
chooses  to  use  an  estimated  average  lot  size,  then  the  user  is  asked  to  enter 
this  value  from  the  terminal.  All  this  new  information  is  passed  back  to  the 
main  program  for  future  use. 

Subroutine  CRDNPC 

This  subroutine  does  intermediate  calculations  for  data  consisting  of  the 
number  rejected  or  the  number  of  nonconformities.  The  number  rejected 
(nonconformities),  the  number  of  data  points,  and  in  the  case  of  data  dealing 
with  nonconformities,  the  lot  size  is  passed  to  this  subroutine.  The  average 
number  rejected  (average  number  of  nonconformities)  is  calculated.  If  the 
data  used  is  the  number  rejected,  then  the  average  fraction  rejected  is 
calculated.  All  this  new  data  is  passed  back  to  the  main  program. 
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Subroutines  CLXP5 ,CLPNPC,CLPC 

These  subroutines  take  the  calculations  passed  to  the  main  program  from 
CRDXRS,  CRDPU,  and  CRDNPC  and  calculate  the  upper  and  lower  control  limits  and 
the  central  line  of  the  control  chart. 

Subroutine  CLXRS 

This  subroutine  calculates  the  control  limits  for  the  X-bar,  R,  and  S 
charts.  The  average  of  X-bar  over  all  the  subgroups,  the  average  range  (std. 
dev.),  the  subgroup  size,  and  the  number  of  subgroups  are  passed  to  this 
subroutine.  The  user  is  given  the  choice  of  having  the  control  limits 
calculated  from  the  data,  using  old  control  limits,  or  using  aimed  at  values. 
If  the  user  chooses  to  have  the  control  limits  calculated  from  data,  then 
these  are  calculated.  If  the  user  chooses  to  use  old  control  limits,  then 
these  are  read  in  from  the  control  limit  data  file,  the  name  of  which  was 
asked  for  at  the  beginning  of  the  main  program.  If  aimed  at  values  for  X-bar 
and  sigma  are  to  be  used,  the  user  is  required  to  enter  these  values  from  the 
terminal  and  then  the  control  limits  are  calculated.  The  control  limits  along 
with  the  data  passed  in  are  passed  back  to  the  main  program. 

Subroutine  CLPNPC 

This  subroutine  is  used  to  calculate  the  control  limits  for  the  np  chart 
or  the  c  chart  or  the  p  chart  if  an  average  lot  size  is  being  used.  The 
values  of  np-bar  and  p-bar  or  c-bar  are  passed  to  this  subroutine.  The  user 
is  given  the  option  of  having  control  limits  calculated  from  the  data,  using 
old  control  limits  or  using  aimed  at  values.  If  aimed  at  values  are  being 
used,  the  user  is  asked  to  enter  them  from  the  terminal.  The  control  limits 
are  then  calculated  the  same  as  if  the  original  data  were  being  used.  If  old 
control  limits  are  being  used,  then  they  are  read  in  from  the  control  limit 
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file*  the  us-*r  alr?oay  set  up.  The  control  limits  are  then  passed  back  to  the 
main  program. 

Subroutine  CLPC 

This  subroutine  calculates  the  control  limits  for  the  u  chart  or  the  p 
chart  (if  using  individual  lot  sizes).  The  values  of  u-bar  or  p-bar  and  the 
number  of  data  points  and  the  number  produced  or  the  number  inspected  are 
passed  to  this  subroutine.  The  user  is  given  the  option  of  having  the  control 
limits  calculated  from  data  or  using  aimed  at  values.  Old  control  limits  are 
not  used  here  since  the  control  limits  depend  on  the  number  produced  or  the 
number  inspected  for  each  data  point.  If  the  user  chooses  to  use  aimed  at 
values,  then  he  or  she  is  asked  to  enter  this  value  from  the  terminal.  Then 
the  control  limits  for  each  pair  of  data  points  are  calculated  and  passed  back 
to  the  main  program. 

Subroutine  Chart 

This  subroutine  plots  the  control  chart  for  X-bar,  R,  s,  p,  np,  and  c. 
The  points  to  be  plotted,  the  control  limits,  and  the  number  of  points  are 
passed  to  this  subroutine.  The  user  is  given  the  option  of  displaying  the 
chart  on  the  screen,  writing  it  to  the  output  file,  or  both. 

Basically,  no  matter  which  chart  is  to  be  printed,  the  procedure  is  the 
same.  The  upper  and  lower  control  limits  and  the  central  line  are  placed  at 
specific  locations  on  the  chart.  The  points  to  be  plotted  are  then  scaled 
between  the  upper  and  lower  control  limits.  If  a  point  falls  above  the  upper 
control  limit  or  below  the  lower  control  limit,  that  point  is  indicated  by  a 
question  mark.  If  it  falls  in  control,  an  asterisk  is  used.  The  values  of 
the  upper  and  lower  control  limits  and  the  central  line  are  printed  on  the 
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chart.  The  user  is  given  the  option  of  saving  the  values  of  the  central  line 

and  control  limits  in  the  control  limit  file  named  in  the  main  program. 

Subroutine  UCHART 

This  subroutine  displays  the  u  or  p  chart  (individual  lot  sizes).  The 
user  is  given  the  option  of  viewing  the  control  chart  on  the  screen  or  having 
it  written  to  the  output  file  or  both.  The  values  for  each  set  of  upper  and 

lower  control  limits,  the  central  line,  and  the  points  to  be  plotted  are 

passed  to  this  subroutine.  This  subroutine  is  basically  the  same  as 
Subroutine  CHART,  except  that  all  values  are  scaled  between  the  highest  valued 
upper  control  limit  and  the  lowest  valued  lower  control  limit.  If  a  point 
falls  outside  its  control  limits,  it  is  indicated  on  the  chart  by  a  question 
mark.  Otherwise,  the  point  is  plotted  as  an  asterisk.  The  user  is  not  given 
the  option  of  saving  these  control  limits  as  they  are  dependent  on  the  lot 
sizes  (number  produced).  If  the  user  wishes  to  use  the  value  of  the  central 
line  again,  he  or  she  can  use  it  as  an  aimed  at  value  in  a  future  run. 

Subroutine  OELXRS,OELPU,DELNPC 

These  subroutines  are  used  to  delete  points  from  the  calculation  of 
control  limits.  These  data  points  are  not  deleted  from  the  permanent  data 
file  or  the  current  data  array  used  by  the  program. 

Subroutine  DELXRS 

This  subroutine  deletes  points  from  the  X-bar,  R,  or  s  charts.  The 
average  of  each  subgroup  and  the  range  (std.  dev.)  of  each  subgroup,  the 
number  of  subgroups,  the  overall  average  and  the  average  range  (std.  dev.)  are 
passed  to  the  program.  The  user  is  asked  to  enter  the  total  number  of 
subgroups  to  be  deleted.  The  sum  total  of  the  averages  of  all  the  subgroups 
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and  Lite  sum  total  of  the  .\.ng,s  (stci.  dev.)  of  all  the  subgroups  arc  thee 
calculated.  The  user  is  asked  to  enter  (one  at  a  time)  the  number  of  the 
subgroup  to  be  deleted.  The  average  and  range  (std.  dev.)  of  the  deleted 
subgroups  are  then  substracted  from  the  sum  totals  and  the  new  overall  average 
and  average  of  the  range  (std.  dev.)  is  calculated.  The  only  values  that 
change  then  are  the  average  of  the  subgroup  averages  and  the  average  of  the 
ranges  (std.  devs.).  Consequently,  when  the  chart  is  printed  out,  all  data 
values  will  be  plotted. 

Subroutine  DELPU 

This  subroutine  deletes  points  from  p  or  u  chart.  The  number  rejected  or 
the  number  of  nonconformities,  the  number  inspected  or  the  number  produced, 
the  average  fraction  rejected  or  the  average  nonconformities  per  unit  and  the 
number  of  data  points  are  passed  to  this  subroutine.  The  user  is  asked  to 
enter  the  total  number  of  pairs  of  data  points  to  be  deleted.  The  total 
number  inspected  or  the  total  number  produced  and  the  total  number  rejected  or 
the  total  number  of  nonconformities  are  calculated.  The  user  is  asked  to 
enter  (one  at  a  time)  the  number  of  the  data  pair  to  be  deleted.  These  values 
are  then  subtracated  from  the  respective  totals  and  the  new  value  of  p-bar  or 
u-bar  is  calculated  and  passed  back  to  the  main  program. 

Subroutine  DELCNP 

This  subroutine  deletes  points  from  the  np  or  c  chart.  The  number 
rejected  or  the  number  of  nonconformities,  the  average  rejected  or  the  average 
number  of  nonconformities  and  in  the  case  of  the  NP  chart,  the  average 
fraction  rejected  and  the  lot  size  are  passed  to  this  subroutine.  The  user  is 
asked  to  enter  the  total  number  of  data  points  to  be  deleted.  The  total 
number  rejected  or  the  total  number  of  nonconformities  is  calculated.  The 
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u:>trr  ib  Lii'.-n  a  Sr.  _u  t.o  cuUi'  th-.  r.u.T.b.r  of  the  dat3  point  to  be  deleted.  The 
values  of  these  data  points  are  then  substracted  from  their  respective  totals 
and  the  new  values  for  the  average  rejected  and  the  average  fraction  rejected 
or  the  average  nonconformities  are  passed  back  to  the  main  program. 

Subroutines  ADDXRS,  ADDPU ,  ADONPC 

These  subroutines  allow  the  user  to  add  new  data  points  to  the  existing 
data  file.  These  points  also  are  added  to  the  permanent  data  file. 

Subroutine  ADDXRS 

This  subroutine  adds  new  data  to  existing  subgroup  data.  The  current 
data  array,  the  number  of  subgroups,  the  subgroup  size,  and  the  name  of  the 
permanent  data  file  are  passed  to  this  subroutine.  The  user  is  asked  to  enter 
the  number  of  subgroups  to  be  added.  If  the  number  of  current  data  points 
plus  the  number  of  new  data  points  is  greater  than  30,  then  enough  subgroups 
are  deleted  from  the  beginning  of  the  file  to  make  room  for  the  new  ones  and 
the  remaining  data  points  are  moved  up  to  the  beginning  of  the  file.  The  user 
is  then  asked  to  enter  each  subgroup  (one  at  a  time).  These  values  are  added 
to  the  current  data  file  and  the  permanent  one.  The  number  or  subgroups  is 
also  updated.  These  new  values  are  passed  back  to  the  main  program. 

Subroutine  ADDPU 

This  subroutine  allows  the  user  to  add  new  data  points  to  the  data  for 
the  P  chart  and  the  U  chart.  The  number  rejected  (number  of  nonconformities) 
and  the  number  inspected  (number  produced)  and  the  number  of  data  pairs  are 
passed  to  this  subroutine.  The  user  is  asked  to  enter  the  number  of  pairs  of 
data  points  to  be  added.  If  the  total  number  of  data  points  (old  and  new)  is 
greated  than  30,  then  the  beginning  data  points  will  be  deleted  to  make  room 


for  tho  new  ones.  Th.  user  will  then  be  asked  to  enter  the  new  data  points 
and  these  will  be  added  to  the  current  data  arrays  and  the  permanent  data 
file. 

Subroutine  ADDNPC 

This  subroutine  allows  the  user  to  add  new  data  points  to  the  data  for 
the  np  and  c  charts.  The  number  rejected  or  the  number  of  nonconformities  and 
the  current  number  of  data  points  are  passed  to  this  subroutine.  The  user  is 
asked  to  enter  the  number  of  data  points  to  be  added.  If  the  total  number  of 
data  points  (old  and  new)  is  greater  than  30,  then  enough  of  the  beginning 
data  points  will  be  deleted  to  make  room  for  the  new  ones.  The  user  will  then 
be  asked  to  enter  the  new  data  points  and  these  will  be  added  to  the  current 
data  array  and  the  permanent  data  file. 

Subroutines  HCDXRS,  HCOPU,  HCDCNP 

These  subroutines  write  to  the  permanent  output  data  file  a  copy  of  the 
raw  data  and  any  intermediate  values  that  are  used  in  plotting  control  charts, 
namely  the  values  of  the  points  that  are  plotted.  The  user  can  obtain  a 
hardcopy  of  this  data  when  the  program  is  completed. 

Subroutine  HISTO 

This  subroutine  plots  a  histogram  of  the  data  used  for  X-bar  and  R  or  s 
charts.  The  raw  data,  the  number  of  subgroups,  and  the  subgroup  size  is 
passed  to  the  subroutine.  The  user  is  then  asked  to  input  the  number  of  cell 
intervals  wanted.  The  number  of  cell  midpoints  and  the  number  in  each  cell 
are  then  calculated  and  the  histogram  is  printed  out.  Each  row  of  asterisks 
in  a  cell  corresponds  to  one  observation  unless  otherwise  indicated  by  the 
scale  on  the  side  of  the  histogram. 
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CONCLUSION 


This  program  is  user  friendly  and  with  little  instruction  can  be  used  by 
anyone  with  a  little  knowledge  in  the  field  of  Quality  Control.  All  of  the 
options  are  designed  to  be  run  independently.  This  causes  some  redundancy 
with  regard  to  the  number  of  subroutines  called,  but  offers  the  user  maximum 
flexibility  when  running  the  program. 

The  Appendix  provides  a  listing  of  the  computer  code. 


UUUUUi 


•••  RAIN  PR OCR AM 

•••  THIS  IS  AN  INTERACTIVE  PR  OCR  ATI  DESIGNED  TO  AID 
•••  IN  THE  DEVELOPMENT  OF  X-BAR.  R.  S.  C.  U.  P,  ANO  NR  CONTROL 
C***  CHARTS 
C 

DIMENSION  DOO.  10).  XBAR  (30) .  RANGEOO I .  SIGMA! 30 ) 
DIMENSION  UOO).  VUCLOO).  VLCLOO).  POO) 

REAL  NIOO).  NROOl.NCOO).  NPOO).  NPB 
INTEGER  SS 
REAL  LCL 

CHARACTER >10  NAME!.  NAMES.  NAMES.  NAMES 

C 

C***  DISPLAY  FIRST  MENU 
C 

URITEO.  1) 

1  FORMAT!  *  1  './//.  10X,  "MENU*.  //> 


URITEO.*)  '  1  ) 
URITEO.*)  *2) 
URITEO.  •)  '3) 
URITEO.  •)  ■  *4) 
READO.  •)  MENU! 


X-8AR  AND  RANGE  CHARTS' 
X-BAR  AND  S  CHARTS' 

C  AND  U  CHARTS* 

P  AND  NP  CHARTS' 


C***  READ  IN  DATA 
C 

URITEO.  •>  'ENTER  THE  NAME  OF  YOUR  DATA  FILE 

1  ((-  10  CHARACTERS)’ 

READO.  3)  NAME! 

3  FORMATIA) 

URITEO.  •>  ‘ENTER  THE  NAME  OF  YOUR  OUTPUT  FILE 

2  ((-  10  CHARACTERS)* 

READO.  3)  NAHE2 

OPEN  (UN  IT*  l*  FILE-NAME!.  STATUS*  'OLD  '  ) 

OPEN (UNIT-2. FILE-NAME2.  STATUS- 'NEW ' ) 

IFIMENUl.EO.  1.0R.MENU1.E0.2)  THEN 

URITEO.*)  'ENTER  THE  NUMBER  OF  SUBGROUPS' 

READO.  *1  NS 

URITEO.*)  'ENTER  THE  SUBGROUP  SIZE' 

READO.  •  I  SS 

READO. •>  ((D(t.  J>.  J-l.SS).  I-l.NS) 

CLOSEIUNIT-l. STATUS* 'KEEP ' ) 

URITEO.*)  'ENTER  THE  NAME  OF  XBAR  CONTROL 
1  LIMIT  FILE«-  10  CHARACTERS)' 

READO,  3)  NAMES 
IFIMENUl.EO.  1)  THEN 

URITEO.  •>  'ENTER  THE  NAME  OF  RANGE  CONTROL 
l  LIMIT  FILE  «-  10  CHARACTERS)’ 

READO.  3)  NAME* 

ELSE 

URITEO.  •)  'ENTER  THE  NAME  OF  SIGMA  CONTROL 
1  LIMIT  FILE  (<  -  10  CHARACTERS)' 

READO.  3)  NAME* 

_ END IF 

EL3EIFIMENU1.E0.4)  THEN 

URITEO.*)  'ENTER  1  IF  WORKING  WITH  FRACTION 
I  REJECTED  <P>' 

URITEO.*)  'ENTER  2  IF  WORKING  WITH  NUMBER 
1  REJECTED  (NP)' 

READO.*)  MNP 

URITEO.*)  'ENTER  THE  •  OF  L0T8  INSPECTED' 

READO.*)  NL 
IF (MNP. EG.  1)  THEN 

URITEO.*)  'ENTER  THE  NAME  OF  CONTROL  LIMIT  FILE' 
READO.  3)  NAMES 
DO  10  I-l.NL 

READ!  1.  •>  NKI).NR(I) 

10  CONTINUE 

ELSE 

URITEO.  •>  'ENTER  THE  LOT  SIZE’ 

READO.  •>  LS 
DO  20  I-l.NL 

READd.*)  NR (  I ) 

20  CONTINUE 

END  IF 

ELSE IF ( MENU 1. EQ. S)  THEN 

URITEO.*)  'ENTER  1  IF  UORKINO  WITH  «  OF 
1  NONCONFORMITIES  (C)’ 

URITEO.  •)  'ENTER  2  IF  UORKINO  UITH  NONCONFORMITIES 
1  PER  UNIT  (U)' 

READO.*)  MCU 

URITEO.  •)  'ENTER  THE  NUMBER  OF  DATA  POINTS' 
READO.*)  N 


non  nnn 


IFIMCU.EO.  1)  IHtN 

URITEO.  •>  'ENTER  THE  NAME  OF  CONTROL  LIMIT  FILE  * 
READ  0.9)  NAMES 
DO  30  IM.N 

READIl.  •)  NCI  I) 

30  CONTINUE 

ELSE 

DO  40  IM.N 

READ!  t.  *)  NCIl).NP(l) 

40  CONTINUE 

END  IF 
END  IF 

49  IFIMENUt.EO.  t.OR.  MENU1.E0.  31  THEN 
•••  DISPLAY  SECOND  MENU 


90 


WRITE  1 9.  90) 
FORMAT!  ‘1  '.  /// 
URITEO.  •)  *t> 
URITEO.*)  *3) 
URITEO.*)  '3) 
IFIMENUt.EO.  I  > 
URITEO.  •) 

ELSE 

URITEO.  •> 
ENDIF 

URITEO.  •)  *9) 
URITEO.*)  *6) 
URITEO.*)  ’7) 
READO.  •)  MENUS 


'  MENU',//) 

HARDCOPY  OF  DATA' 
HISTOGRAM' 

X-BAR  CHART' 

THEN 

'4)  R  CHART' 

‘4)  S  CHART' 

ADD  NEW  SUBCROUPS' 
DELETE  SUBCROUPS' 
EXIT  PROGRAM' 


•••  CALL  APPROPRIATE  SUBROUTINE  ACCORDING  TO  OPTION  CHOSEN 

IFIMENU2.E0.  I.  AND.  MENU  I.  EO.  II  THEN 

CALL  CRDXRSID.  NS.  9S.  XBAR#  RANGE.  XBB.  RANB.  MENU l  > 

CALL  HCDXRSID.  NS.  SS.  XBAR.  RANGE.  HENU1.  *49) 
ELSEIFIMENU2.EQ.  1 .  AND.  HENU1 .  EQ.  2 >  THEN 

CALL  CRDXRSID.  NS,  SS.  XBAR.  SIGMA,  XBB.  SIGB.  MENUI  > 

CALL  HCDXRSID.  NS.  SS.  XBAR,  SIGMA.  MENUI.  *49> 

ELSE IF | MENUS. EG. 2)  THEN 

CALL  HISTOID.  NS.  SS,  *49) 

ELSEIFIMENU2.ED. 3. AND. MENU!. EQ. U  THEN 

CALL  CRDXRSID.  NS.  SS.  XBAR.  RANGE.  XBB.  RANB.  MENUI  ) 

CALL  CLXRSIMENU1,  MENUS.  XBB.  RANB.  NS.  SS.  UCL.  LCL. 

CL, NAMES) 

CALL  CHART  I  XBAR.  UCL.  LCL.  CL.  NS.  NAMES.  MENUI .  MENUS. 
MNP,  MCU.  *49) 

ELSEIFIMENU3.E0.  3.  ANO.MENUl.Ea.  2)  THEN 

CALL  CRDXRSID.  NS.  SS. XBAR.  SIGMA. XBB. SIGB.  MENUI ) 

CALL  CLXR8IMENU1.  MENUS.  XBB.  SIGB.  NS.  PS.  UCL.  LCL. 

CL.  NAMES) 

CALL  CHART  I  XBAR.  UCL.  LCL.  CL.  NS.  NAME3.  MENUI.  MENUS. 
MNP. MCU.  *49) 

ELSEIFIMENUa.  EG.  4.  AND.  MENUI .  EO.  1  >  THEN 

CALL  CRDXRSID. NS.  SS. XBAR. RANGE. XBB. RANB.  MENUI ) 

CALL  CLXRSIMENU1.  MENUS.  XBB.  RANB.  NS.  SS.  UCL.  LCL. 
CL.NAME4) 

CALL  CHART IRANGE.  UCL.  LCL.  CL.  NS.  NAME4.  MENUI.  MENUS. 
MNP. MCU.  *49) 

ELSE IF I MENUS. EQ. 4. AND. MENUI. EG. 2)  THEN 

CALL  CRDXRSID.  NS.  SS.  XBAR.  SIGMA.  XBB.  SIGB.  MENUI  I 
CALL  CLXR8! MENUI.  MENU2.  XBB.  SIGB.  NS.  SS.  UCL.  LCL. 

CL.  NAME4  ) 

CALL  CHART  I  SIGMA.  UCL.  LCL.  CL.  NS.  NAME4.  MENUI.  MENU2. 
MNP.  MCU.  *49) 

EL8EIF IMENU2. EO. 9 )  THEN 

CALL  ADDXRSI D,  NS.  8S.  NAME l .  *49) 

ELSE IF (MENUS. EQ. &)  THEN 

URITEO.*)  'ARE  THESE  SUBGROUPS  TO  BE  DELETED  FROM' 
URITEO.*)  '1)  X-BAR  CHART' 

IFIMENUI.EC. I )  THEN 

URITEO.*)  '2)  R  CHART' 

ELSE 

URITEO.*)  '2)  S  CHART' 

ENDIF 

READO.*)  MDEL 
IFIMDEL.EO.il  MENU2-3 
IFIMDEL. EO. 2)  MENU2-4 
IF  I  MDEL.  EO.  I.  AND.  MENUI.  EO.  I )  THEN 

CALL  CRDXRSID.  NS.  SS.  XBAR.  RANGE.  XBB.  RANB.  MENUI ) 
CALL  OELXRSIXBAR.  RANGE.  NS,  XBB.  RANB.  MENUI) 

CALL  CLXRSIMENU1.  MENUS.  XBB.  RANB.  NS.  SS.  UCL.  LCL. 
CL.  NAMES  > 


) 


or>r> 


K-:v. 


CALL  CHART <  XBAR.  UCL,  LCL.  CL.  NS.  NAHE3.  HENUl  .  HENU2. 

I  HNP.  HCU.  *45) 

ELSE  IF  <  HOEL. EQ. 1 . AND. HENUl . EQ.  2)  THEN 

CALL  CRDXRS ( D. NS. SS.  X8AR.  SIGHA.  XBB.  SIGB.  HENUl  ) 
CALL  DEL  XRS(XBAR,  SICPIA.  NS.  XBB.  SI  GB.  MENU  I  > 

CALL  CLXRS(flENUl.  HENU2.  XBB.  SICB.NS.  SS. UCL.  LCL. 

1  CL.NAHE3) 

CALL  CHART ( XBAR. UCL. LCL. CL.  NS. NAME 3. MENU l .  HENU2. 

1  O.  O.  *45) 

ELSEIFIHOEL.EO. 2. AND.HENU1.EQ. 1 >  THEN 

CALL  CRDXRSCD. NS. SS.  XBAR. RANGE. XBB. RAN8.  MENU! > 
CALL  DELXRSt  XBAR. RANGE. NS. XBB. RAN8.  HENUl ) 

CALL  CLXRSIHENUl .  HCNU2.  XBB.  RANB.  NS.  SS.  UCL.  LCL. 

I  CL. NAME4 I 

CALL  CHART  ( RANGE.  UCL.  LCL.  CL.  NS.  NAM&4.  HENUl .  HENU2. 
■  O.O.  *43) 

ELSEIF<HDEL. EQ. 2.  AND. HENUl. EQ. 2)  THEN 

CALL  CRDXRSCD.  NS.  SS. XBAR. SIGHA. XBB.  SIGB.  HENUl ) 
CALL  DELXRSl XBAR, SIGHA. NS. XBB. SIGB.  HENUl ) 

CALL  CLXRS( HENUl »  HENU2.  XBB.  SICB.NS.  SS.  UCL.  LCL. 

1  CL. NAHE4) 

CALL  CHART (SIGHA. UCL. LCL. CL. NS.  NAHE4. HENUl . MENU 2, 
1  0.0.  *43) 

END  IF 

ELSEIF(HENU2.EQ. 7)  THEN 
GO  TO  lOOO 
END  IF 

***  DISPLAY  SECOND  HENU 

ELSE IF « HENUl.  EQ. 3. OR.  HENUl .EQ. 4)  THEN 
WRITEO.  60) 

60  FORHATC  *1 '•///•  '  HENU'./) 

WRITE <3. •)  '1)  HARDCOPY  OF  DATA' 

IF (HENUl. EO. 4)  THEN 

IF( HNP. EQ. 1 >  THEN 

WRXTE(3. *>  *2)  P  CHART' 

ELSE 

WRITE O.  •  )  '2)  NP  CHART' 

END  IF 

ELSE 

IF (HCU. EQ. 1 )  THEN 

WRITEO.  •)  *2)  C  CHART' 

ELSE 

WRITE <3.  *)  '2)  U  CHART' 

END  IF 


ENDIF 
WRITE (3.  *) 
WRITEO.  *) 
WRITEO.  *) 


'3>  ADD  HORE  DATA' 
'4)  DELETE  DATA' 
'3>  EXIT  PROGRAH' 


READO.  •  >  HENU2 

C 

C***  CALL  APPROPRIATE  SUBROUTINE  ACCORDING  TO  OPTION  CHOSEN 

I F ( MENU2. EQ. 1 )  THEN 

IF  <  HENUl • EQ. 3. AND. HCU. EQ. I )  THEN 
CALL  HCDCNPCNC,  N.  0.  *43) 

ELSEIF (HENUl .  EQ. 3. AND*  HCU. EQ. 2)  THEN 

CALL  CRDPU(NC«  NP.  N#  HENUl.  U.  UB.  0*  0.  HENU2) 

CALL  HCDPU<NC.  NP.  U.  N.  HENUl.  *43) 

ELSE IF ( HENUl. EQ, 4.  AND. HNP.  EQ.  1 )  THEN 

CALL  CRDPU(NR.Nt.NL*HENUI.P.  PB.HENP.  ALS.HCNU2) 
CALL  HCDPUCNR.NI.P.NL.  HENUl.  *43) 

ELSEIF ( HENUl .  EQ.  4.  AND.  HNP.  EQ.  2)  THEN 
CALL  HC0CNP(NR.  NL.LS.  *43) 

ENOtF 

ELSE 1F(H£NU2. EQ. 2. AND. HENUl. EQ. 3)  THEN 
1FCHCU.EQ. I)  THEN 

CALL  CRDNPCINC.N.  O.ANC.O) 

CALL  CLPNPC (HENUl.  O. ANC. 0.  UCL. LCL. CL.  NAHE3) 
CALL  CHART (NC.  UCL.  LCL.  CL.  N.  NAME 3.  HENUl. 

I  HENU 2. O. HCU.  *43) 

ELSE 

CALL  CRDPU(NC.NP.N.  HENUl.  U.UO.  0.0.  HENU2) 

CALL  CLPU(HENUi .  N.  NP.  UB.  VUCL.  VLCL.  CL  > 

CALL  UCHARTCU.  VUCL. VLCL.  CL.  N.  HENUl . *43) 

ENDIF 

ELSEIFCHCNU2. EQ. 2. AND. HENUl. EQ. 4 )  THEN 
IFIMNP.EQ.  1>  THEN 

CALL  CRDPUXNR.Nl.  NL.  HENUl.  P.  PB.  HENP.  ALS.HENU2) 
IF(HENP.EQ. 1)  THEN 

CALL  CLPU< HENU1.NL.  NI, PD.  VUCL.  VLCL.  CL) 

CALL  UCHARTfP.  VUCL.  VLCL.  CL. NL.  HENUl.  *43) 


v 
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CALL  CLPNPC ( MENU l . MNP. PB.  ALS.  UCL.  LCL.  CL.  NAMES ) 
CALL  CMART(P. UCL. LCL. CL. NL.  NAMES.  MENUI. 

MENUS, MNP. O. *45) 

END  IF 

ELSE 

CALL  CRDNPC (NR.  NL.  LS.  NPB.  PD ) 

CALL  CLPNPC < MENU 1, MNP. PB. NPB.  UCL.  LCL.  CL.  NAMES > 

CALL  CHART < NR, UCL. LCL. CL. NL. NAMES.  MENU1. 

1  MENU2. MNP. O. *45) 

END  IF 

ELSE I F ( MENU2. EQ. 3. AND. MENU 1 .  EG.  3  >  THEN 
IF(MCU. EQ. 1)  THEN 

CALL  ADDCNPINC. N. NAME I .  *45) 

ELSE 

CALL  ADDPUINC. NP, N  MENU l.  NAME 1. *45) 

END  IF 

ELSE IF ( MENU2. EQ. 3. AND. ME'  EO.  4 >  THEN 

IF ( MNP. EO. I )  THEN 

CALL  ADDPUINR. NI. NL. MENU1.  NAME1.  *45) 

ELSE 

CALL  ADDCNPINR.  NL.  NAME I .  *45) 

END  IF 

ELSE IF ( MENU2. EQ. 4. AND.  MENU I . EO.  3 )  THEN 
IF (MCU. EQ. 1 )  THEN 

CALL  CRDNPC CMC.  N. O. ANC, 0) 

CALL  DELCNPCNC.  ANC. O. N. 01 

CALL  CLPNPC ( MENU 1 . 0. ANC. 0. UCL. LCL.  CL.  NAMES) 

CALL  CHARTCNC. UCL. LCL. CL. N. NAMES.  MENU!. 

!  MENUS.  O.  MCU.  *45) 

ELSE 

CALL  CRDPU(NC.  NP.  N.  MENU1.  U.  UD.  O.  O*  MENU2) 

CALL  DELPUINC.  ANC.NP.N) 

CALL  CLPU < MENUI.  N.  NP.  U&.  VUCL.  VLCL.  CL > 

CALL  UCMARTIU. VUCL. VLCL. CL.  N.  MENU1.  *45) 

END  IF 

ELSEXFi MENUS. EQ. 4. AND. MENU1. EQ, 4)  THEN 
IF (MNP.  EQ. 1 )  THEN 

CALL  CRDPU (NR.  NI , NL. MENU I.  P.  PB.  MENP.  ALS.  MENU2) 

CALL  D£LPU(NR.PB. NI.  NL) 

IF C  MENP. EO.  1 )  THEN 

CALL  CLPU ( MENU  1.  NL.  N I.  PB.  VUCL.  VLCL.  CD 
CALL  UCHARTfP. VUCL.VLCL.CL.NL.  MENUl.  *45) 

ELSE 

CALL  CLPNPC  < MENU l . MNP.  PB.  ALS.  UCL.  LCL.  CL > 

CALL  CHART (P.  UCL. LCL.  CL.  NL.  NAMES.  MENU1. 

1  MENUS.  MNP.  O. *43) 

END  IF 

ELSE 

CALL  CRDNPC! NR. NL.  LS.  NPB.  PB ) 

CALL  DELCNP (NR. NPB. PB. NL.  LS) 

CALL  C LPNPCC MENU 1. MNP. PB.  NPB.  UCL.  LCL.  CL.  NAMES) 

CALL  CHART (NR. UCL. LCL. CL. NL. NAME3.  MENU!. 

I  MENUS.  MNP.  O.  *45) 

END  IF 

ELSE  IF!  MENUS.  EQ.  3)  THEN 
00^1000 
ENDIF 
END  IF 
lOOO  STOP 
END 

•••  THIS  SUBROUTINE  DELETES  DATA  POINTS  FROM  THE  X— BAR.  R. AND  S  CHARTS 

SUBROUTINE  DEL XRS( A.  B» NS.  AB.  BB.  MENU 1 ) 

DIMENSION  A (SO) .  8(30) 

WRITE! 3.*)  'ENTER  THE  TOTAL  NUMBER  OF  SUBGROUPS  TO  BE 
1  DELETED ' 

READO.  •)  NO 
NNS-NS-ND 
BB-BB*NS 
A8-AB*NS 
DO  10  1*1.  ND 

WRITE (9.*)  'ENTER  SUBGROUP  NUMBER  TO  BE  DELETED' 

READ! 3. •)  K 

AB-AB-AOO 

BB-BB-B(K) 

IO  CONTINUE 
AB-AB/NNS 
B8-BB/NNS 
RETURN 


uuu 


C***  THIS  SUBROUTING  DELETES  DATA  POINTS  FROM  THE  P  AND  U  CONTROL  CHARTS 
C 

SUBROUTINE  DCLPUC A.  AD. D.  N > 

DIMENSION  A<  30 ) .  8  C 30 ) 

WRITE C3.  *>  'ENTER  THE  TOTAL  •  OF  PAIRS  OF  DATA  POINTS 
1  TO  BE  DELETED* 

READ <3. *)  ND 
TLS-O.  0 
DO  3  1*1.  N 

TLS—TLS+B ( I > 

3  CONTINUE 
A8-AB*TLS 
DO  10  1*1 »  ND 

WRITE <3. *)  'ENTER  DATA  POINT  •  TO  BE  ELIMINATED 
READC 3.  •)  K 

AO-AB-A(K) 

TLS-TLS— B IK ) 

10  CONTINUE 
AB-AB/TLS 
RETURN 
END 

***  THIS  SUBROUTINE  DELETES  POINTS  FROM  THE  C  AND  NP  CONTROL  CHARTS 

SUBROUTINE  DELCNPCA. AB. BB*  N. LS) 

DIMENSION  AC30) 

WRITE (3. *1  'ENTER  THE  TOTAL  #  OF  DATA  POINTS  TO  BE  DELETED' 

READC 3* • >  ND 
NNS-N-ND 
AB-AB*N 
DO  10  1*1* N 

WRITE (3.*)  'ENTER  DATA  POINT  •  TO  BE  DELETED' 

REA0C3.  •  )  K 
AB—AS— ACK ) 
lO  CONTINUE 

BB*AB/ ( NNS*LS ) 

AB-A8/NNS 

RETURN 

END 

***  THIS  SUBROUTINE  PRINTS  A  HISTOGRAM  TO  THE  SCREEN 
***  AND  TO  AN  OUTPUT  FILE 

SUBROUTINE  HISTO CD.  NS. SS.  •) 

DIMENSION  DC 30.  10).  KC0UNTC20).  NC0UNTC20).  NCOORC 12) 

DIMENSION  CMC 20) 

INTEGER  SS.C2 
CHARACTER *63  COL. DASH 
CHARACTER*70  LONG 
WRITE!  3.*)  '  MENU' 

WRITE! 3. *)  '  ' 

WRITE! 3.*)  'D  PRINT  HISTOORAM  TO  SCREEN' 

WRITE!3. *)  '2)  PRINT  HISTOORAM  TO  FILE' 

WRITE! 3.*)  '3)  BOTH' 

READ! 3.  *)  MENUH 
DMA  X -DU.  1) 

DMIN-DC 1.  I) 

DO  20  1*1.  NS 

DO  10  U-l.SS 

IF C DC  I,  J).LT.  DMIN)  DMIN-DCI.J) 

IF  (DC  I,  J).  CT.  DMAX  )  DMAX-DCI.J) 

10  CONTINUE 

20  CONTINUE 

DO  23  1-1.63 

COLCI:  I)*'*' 

DASH C I  I)*'-' 

23  CONTINUE 

WRITE! 3. • )  '  ' 

WRITE! 3.*)  'INPUT  •  OF  CELL  INTERVALS  CC-20)' 

READ!  3.*)  C2 
CI*(  DMAX -DMIN )/C2 
DO  30  J-l. C2 

KCOUNTC J)— O. O 
30  CONTINUE 

DO  40  1-1.  NS 

DO  33  K-l.SS 

J-INT!  CD!  I.K)-DMIN)/CI  >«-1.0 
IFCU.GT.20)  U-20 
KCOUNTC JJ-KCOUNT C  J) *1 . O 
33  CONTINUE 

40  CONTINUE 
KMAX-O 
DO  43  I-1.C2 


M\ 

L- 


1 F  ( KCOUNT  ( I ) .  GT.  KMAX  >  UMAX -KCOUNT ! I ) 

43  CONTINUE 

IF!M£NUH.  EQ. 2)  GQ  TO  90 
NDUM-KHAX/20M 
FACTOR-NDUM 
DO  30  I-I. 20 

NCOUNT ( I> - ( KCOUNT  < I ) ♦. 3  > /FACTOR 

:>o  continue 

WRITE (3. 33) 

33  FORMAT  i  '  l  ' ,  20X.  'HISTOGRAM  *  ) 

DO  80  1-20.  I.  -| 

LONG ! |:  !)■'•* 

DO  60  K-2. 70 

LONG ! K:  K ) «  '  * 

60  CONTINUE 

IPT-3 

DO  63  J-l.  20 

IF< I.CT. NCOUNT <J>>  THEN 
IPT-IPTO 

ELSE 

LONG! 1PT :  IPT>«'*' 

TPT-IPTM 

LONG< IPT :  IPT)-'*' 

IPT-JPT+2 
END  IF 

65  CONTINUE 

IID-FACTQR*I 
lF(nOD< I. 3) . EO.  O)  THEN 
WRITE<3.70)  1 1 D. LONG 
70  FORMAT !  *  '.13.  A70) 

ELSE 

WRITE! 3* 73)  LONG 
73  FORMAT !  *  '.3X.A70) 

END  IF 

80  CONTINUE 

WRITE! 3. 83)  DASH!1:63) 

83  FORMAT!'  '.3X.A6S) 

WRITEO.  *)  'PRESS  1  TO  CONTINUE' 

READ! 3#  * )  ND 
IFCMENUH. EQ. I)  GO  TO  130 
C 

C***  PRINT  HISTOGRAM  TO  OUTPUT  FILE 

C 

90  WRITE! 2.  93 > 

93  FORMAT!  '1  30X.  'HISTOGRAM'. /> 

WRITE! 2. 10OI 

100  FORMAT! 'O'.  30X.  'FREQUENCY') 

WRITE! 2. 103) 

103  FORMAT!'  'CELL  MXDPT.  ') 

DASH!!: !>-'♦' 

DASH! 63:  63>-'«-' 

WRITE! 2.  1 10)  DASH! I: 63) 

HO  FORMAT!  12X.63A) 

NDUM-KMAX/60+1 
FACTOR -NDUM 
DO  133  I— I . C2 

K-IKCOUNT! I )♦. 31/FACTOR 

CM! I )-DMlN+<2.  O*! I-I >*1 .  O ) *C I /2. O 

IF!K. EQ. 0)  THEN 

WR!TE!2. 113)  CM! 1 1 

113  FORMAT!  12X.  '!  '/'  F10.  2.  1  X.  ' !  V12X.  ' !  '  > 

ELSE 

WRITEI2, 120)  COL ! 1 : K ) 

120  FORMAT! 12X* A) 

WRITE!2.  123)  CM! I ) .  COL ! 1: K ) 

123  FORMAT!'  '.  F10.  2.  IX,  A) 

WRITE! 2. 120)  COL ! 1 : K ) 

WRITE12. 130) 

130  FORMAT!  12X.  '♦') 

END  IF 

133  CONTINUE 

WRITE!2.  120) 

WRITE!2.  120)  DASH!!:  63) 

DO  140  l-l. 12 

NCOORi X )-I*FACT0R»3 
140  CONTINUE 

WRITEI2.  143)  !NC00R!X),  I-l,  12) 

143  FORMAT! 12X.  1213) 

130  RETURN  1 
ENO 

r* 

THIS  SUBROUTINE  CALCULATES  X-BAR.  THE  RANGE  OR  THE 
STANDARD  DEVIATION  FOR  EACH  SUBGROUP  AND  ALSO 


-.'S', 


,V 


J  J 


X-DOUBLE  BAR.  R-BAR.  AND  S-BAR 

SUBROUTINE  CRDXRS <  D.  NS.  SS.  AI.A2.  AID.  A2B. MENU1 > 

DIMENSION  DOO.  10).  Al  (30) .  A2O0) 

INTEGER  SS 
AIB-O. O 
A2B-0.0 

IF (MENU! . EQ. 1 )  THEN 
DO  20  I * 1 . NS 
A I  (  I  )-0.0 
XHI-D( I.  I  ) 

XLO*D( I.  1  ) 

DO  10  J* l . SS 

Al ( I )«Al ( I >♦□( I*  J) 

IF<D( I.  J).LT. XLO)  XLO-D( I . J) 

IF ( D  <  I. J ) . CT.  XHI )  XMl-D(I.J) 

CONTINUE 

Al ( I )-Al < I )/SS 

A1B-A1B+A1 ( I ) 

A2 ( I ) “XHI -XLO 
A2B-A2B*A2( I > 

CONTINUE 

A10-A1B/NS 

A2B-A2B/NS 

ELSE 

DO  30  1*1. NS 
AKD-0.0 
DO  30  J*  l .  SS 

Al ( I )*A1 ( I )*0( I* J) 

CONTINUE 

ai<i  >-ai(  n/ss 

A 1 B»A 1 B ♦A 1 ( I ) 

A2< I )*0. O 
DO  40  J-I.  SS 

A2( I )*<  D( I.  J)-A1 (I > >**2*A2< I ) 

CONTINUE 

A2( I )*S0RT(A2< I >/<SS-l ) ) 

A2B-A2B ♦A2 ( 1 ) 

CONTINUE 
Al B*AIB/NS 
A2B*A2B/NS 
END  IF 
RETURN 
END 

THIS  SUBROUTINE  CALCULATES  THE  FRACTION  REJECTED  OR 
NONCONFORMITIES  PER  UNIT  FOR  EACH  DATA  PAIR  AND 
THE  AVERAGE  FRACTION  REJECTED  AND  THE  AVERAGE 
NONCONFORMITY  PER  UNIT 

SUBROUTINE  CRDPUCCDl*  CD2.  N. MENU! *  CDB. CDB8. MENP.  ALS.  MENU2) 

DIMENSION  C0K30).  CD2O0).  CDB (30) 

CDBB-0.0 
CR-O.  O 
CN-O.  O 
DO  20  1*1.  N 

CDB (I >-0.0 

CDBC I )*CDi < I ) /C02< I ) 

CR-CR+CDMI) 

CN-CN*CD2< I > 

CONTINUE 

CDBB-CR/CN 

IF (MENU2. EQ.  1 )  GO  TO  40 
IF(MENUl.EO.  4)  THEN 

WR1TEC3.  •>  '  DO  YOU  WISH  TO  USE ' 

WR ITE<  3#  *  I  'll  INDIVIDUAL  LOT  SIZES' 

WRITE (3.*)  '2>  ESTIMATED  AVERAGE  LOT  SIZE' 

WRITE (3. •>  #3>  CALCULATED  AVERAGE  LOT  SIZE' 

READ! 3.  *)  MENP 
IF  (MENP.  EQ.  2)  TtCN 

WRITE! 3. •>  'ENTER  ESTIMATED  AVERAGE  LOT  SIZE' 

READ (3. »>  ALS 
ELSE IF (MENP. EO. 3)  THEN 
DO  30  I-l.N 

ALS— AL3+CD2 ( I > 

CONTINUE 
ALS-ALS/N 
END  IF 
END  IF 
RETURN 
END 

THIS  SUBROUTINE  CALCULATES  THE  AVERAGE  NUMBER  REJECTED  AND  AVERAGE 


FRACTION  REJECTED  OR  THE  AVERAGE  NUMBER  OF  NONCONFORMITIES. 

SUBROUTINE  CRDNPC CCRD. N. LS. CNB. CPO ) 

DIMENSION  CROC 30 ) 

CNQ-0.0 
DO  10  I-l.N 

CNB— CNB ♦CRDC I ) 

CONTINUE 

CNB-CNB/N 

IF(LS.NE.O)  CPB-CNB/LS 

RETURN 

END 

THIS  SUBROUTINE  ALLOWS  THE  USER  TO  ADD  NEW  SUBGROUPS  TO 
THE  EXISTING  DATA  FILE 

SUBROUTINE  ADDXRSCD.  NS.  SS.  NFILE.  « > 

DIMENSION  DC 30.  10) 

INTEGER  SS 
CHARACTER*10  NFILE 

WR I  TEC  3.* I  'ENTER  THE  NUMBER  OF  SUBGROUPS  TO  BE 
1  ADDED* 

READ (3. •  >  NA 
NT— NA*NS 

IFCNT.CT. 30)  THEN 
ND-NT-30 
DO  20  I-l. NS-ND 
DO  10  J-l.SS 

DC  I.  JJ-DCND+I.  J) 

CONTINUE 

CONTINUE 

NN-N5-ND+1 

NS-30 

ELSE 

NN— NS+1 
NS-NS^NA 
END  IF 

DO  40  I«NN,  NS 
K-I-CNN-1 ) 

WRITEC3.  30)  K 

FORMAT ( '  'ENTER  SUBGROUP  ft  '.  12.  '  TO  BE  ADDED') 
READC3.  #)  CDC  I.  J).  J-l.  SS) 

CONTINUE 

CLOSE C UNIT-1 ,  STATUS- 'DELETE ' ) 

OPEN C UN I T— 1 . NAME-NFILE.  STATUS- 'NEW ' ) 

WRITEC1.  #>  C  CDC  I.  J).  J-l,  SS).  I-l.NS) 

RETURN  I 
END 

THIS  SUBROUTINE  ALLOWS  THE  USER  TO  ADD  DATA  POINTS 
TO  THE  EXISTING  DATA  FOR  A  C  OR  NP  CONTROL  CHART 

SUBROUTINE  ADDCNP C A. N.  NFILE. *> 

DIMENSION  AC30) 

CHARACTER* 10  NFILE 

WRITEC  3. *)  'ENTER  THE  ft  OF  DATA  POINTS  TO  BE  ADDED' 

READC3.  *)  NA 
NT-NA+N 

IFCNT. OT. 30)  THEN 
ND-NT-30 
DO  10  1-1. N-ND 
ACD-ACND+1) 

CONTINUE 
NN-N-ND* 1 
N-30 

ELSE 

NN— N*I 
N-N+NA 
ENDIF 

DO  30  I-NN. N 
P- I -CNN-1 ) 

WRITE (3. 20)  K  _  _ 

FORMAT <  '  'ENTER  DATA  POINT  ft  ',12.  '  TO  BE  ADDED') 
READC  3.  * )  AC  I) 

CONTINUE 

CLOSECUNIT-l,  STATUS- 'DELETE ' > 

OPEN (UNIT- 1.  NAME-NFILE.  STATUS- 'NEW ' ) 

DO  40  I-l.N 

WRITEC  1 .  *  )  Ad) 

CONTINUE 
RETURN  1 
END 
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THIS  SUBROUTINE  ALLOWS  THE  USER  TO  ADD  DATA  POINTS 
TO  THE  EXISTING  DATA  FOR  A  P  OR  U  CHART 

SUBROUTINE  AODPUiA.  B.  N.  MENU! .  NF1LE.  * ) 

D I HENS I ON  A!  30). 8(30) 

CHARACTER* 10  NFILE 

WRITE! 3. *)  *  ENTER  THE  4  OF  PAIRS  OF  DATA  POINTS  TO  BE 
1  ADDED ' 

READ! 3.  *)  NA 
NT-NA+N 

IF! NT.  CT.30)  THEN 
ND-NT-30 
DO  10  1*1  *  N—ND 
A! I ) -A!ND+I ) 

B! I )*B (ND*1 ) 

CONTINUE 

NN-N-ND*! 

N-30 

ELSE 

NN-N+1 
N-N+NA 
END  IF 

DO  40  I-NN. N 

IF ( MENU 1.  EQ. 3)  THEN 
WRITEI3. 20) 

FORMAT!'  'ENTER  •  OF  NONCONFORHI TIES  AND  UNITS 
I  PRODUCED') 

ELSE 

WRITE!  3.  30) 

FORMAT!'  '.  'ENTER  •  OF  REJECTS  *  •  INSPECTED*) 
END  IF 

READ !  3.  *)  A!  I  ).  B  <  1  ) 

CONTINUE 

CLOSE  <  UN I T« 1 . STATUS- ' DELETE ' > 

OPEN ( UN I T-l.NAME-NF I LE.  6TATUS- 'NEW ' ) 

IF ! MENU 1 • £0.  3)  THEN 
DO  30  t-l. N 

WRITE! 1.  *)  A!I).B!I> 

CONTINUE 

ELSE 

DO  60  I-l.N 

WRITE!!.*)  B!I).A(I) 

CONTINUE 
ENDIF 
RETURN  1 
END 

THIS  SUBROUTINE  WRITES  A  HARDCOPY  OF  THE  RAW  DATA 
AND  THE  POINTS  TO  BE  PLOTTED  FOR  AN  X-BAR.  R.  OR  S 
CONTROL  CHART 

SUBROUTINE  HCDXRS!D«  NS#  SS.  A*  B.  MENU1. • ) 

DIMENSION  D!30«  10).  A!30>.  B ! 30) 

INTEGER  SS 
WRITE! 2. 10) 

FORMAT! '1 './//.  33X.  'SU0OROUP  DATA'./,  'O'.  IX.  'S0N'.33X. 

1  'RAW  DATA'.  //) 

DO  20  I-l.NS 

WRITE(2.  13)  I,  !D!X. J).  J-l, SS) 

FORMAT!'  '.  IX.  12.  2X.  10!  IX.  F6.  2)) 

CONTINUE 

IF!MENU1.E0.  1)  THEN 
WRITE! 2* 30) 

FORMAT! '1 ',///.  33X,  'SUBGROUP  DATA './//,  33X.  'XBAR  '. 

1  9X.  'RANGE',  //) 

00  30  I-l.NS 

WRITE <2. 40)  A! I ) . B< I ) 

FORMAT!'  '.  28  X,  F9.  4,  3X.  F9.  4) 

CONTINUE 

ELSE 

WRITE! 2.  60) 

FORMAT!  '1  ',///.  33X,  'SUBGROUP  DATA'. ///.  33X,  'XBAR'. 

1  MX,  'S',  //) 

DO  80  1-1, NS 

WRITE(2.  70>  A ! I ) »  B (  I  ) 

FORMAT!'  '.  28X.F9.  4.  3X.F9.4) 

CONTINUE 
ENDIF 
RETURN  1 
END 

THIS  SUBROUTINE  WRITES  A  HARDCOPY  OF  THE  RAW  DATA  AND 
THE  POINTS  TO  BE  PLOTTED  FOR  A  C  OR  A  HP  CONTROL  CHART 
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SUBROUTINE  HCDCNPINHD.  NOP.  LS.  •> 

REAL  NH0I3OI 
IFILS.EO.  0>  THEN 
WRITEI2. 10) 

0  FORMAT <  '1  ///.  32 X.  'DATA  SET  FOR  C -CHART  '  > 

WRITE  12.  30) 

D  FORMAT  1///3BX.  'NUMBER'.  /40X,  'OF'.  /34K.  'NONCONFORMITIES' 

I  //) 

DO  30  1*1.  NOP 

WRITE!  2.  23)  NHD(I) 

»  F0RMATI37X.F7. 2) 

3  CONTINUE 

ELSE 

WRITE (2*  40) 

FORMAT!  '1  ',///.  32X.  'DATA  SET  FOR  NP-CHART') 

WRITE  1 2.  49) 

FORMAT! 'O'.  34X.  'LOT ',  8X.  'NUMBER'/'  '.  33X,  'SI2E'. 

1  7X,  'REJECTED'. //) 

DO  99  I-l.NDP 

WRITE  (2.  90)  LS.NH01I) 

FORMAT!'  '.  33X.  14.  7X.  F7.  2) 

CONTINUE 
ENOIF 
RETURN  I 
END 

THIS  SUBROUTINE  WRITES  A  HARDCOPY  OF  THE  RAW  DATA  AND 
THE  POINTS  TO  BE  PLOTTED  FOR  A  P  AND  U  CHART. 

SUBROUTINE  HCDPUIHHD.  HH02.  HH01. NDP. MENU1.  •) 

DIMENSION  M1O!3O>.HHO2!3>.KHD1!30> 

IF! MENU 1 • EG. 3 )  THEN 
WRITE  (2.  10) 

FORMAT! '1 '. ///32X.  'DATA  SET  FOR  U-CHART ' ) 

WRITE <2.  20) 

FORMAT!///.  19X.  'NUMBER'.  14X.  'NUMBER*.  8X. 

1  'NONCONFORMITIES'/'  '.  20X.  'OF '.  19X.  'OF  UNITS '. 

2  13X.  'PER'/'  '.  14X.  'NONCONFORMITIES'.  BX.  'PRODUCED'. 

3  12X.  'UNIT'.  //) 

DO  40  1-1.  NDP 

WRITE <2.  30)  HHD2!  I  >.  HHD!  I  ).  HHD1  1 1  > 

FORMAT!'  '.  16X.F7.2.  13X.F7.2.  10X.F9.4) 

CONTINUE 

ELSE 

WRITE12.  90) 

FORMAT!  '1 '.  ///.  32X.  'DATA  SET  FOR  P-CHART') 

WRITE12.  60) 

FORMAT!///.  29X.  'NUMBER'.  IOX.  'NUMBER'.  IOX.  'FRACTION', 

1  /'  '.22X.  'INSPECTED'.  7X.  'REJECTED'.  IOX.  'REJECTED'.//) 

DO  BO  I-l.NDP 

WRITEI2,  70)  HHD2!  I ).  HHD!  1 ), HMD!  !  I  > 

FORMAT!  '  ',  23X.  F7.  2.  SX.  F7.  2.  9X.  F9.  4) 

CONTINUE 
ENDIF 
RETURN  I 
END 

THIS  SUBROUTINE  CALCULATES  THE  CONTROL  LIMITS  FOR  A  X-BAR. 

R.  OR  S  CONTROL  CHART 

SUBROUTINE  CLXRS!MENUI.  MENU2.  XBB.  SB.  NS.  SS.  UCL.  LCL.  CL.  NFILE) 
DIMENSION  A2I10).A1!IO).D4!10).D3I10).  84!  10).  B3!  10) 
DIMENSION  D2!10).D1!  10)  .SD2!  10).B2!10).B1!10).C2!10) 

INTEGER  SS 
REAL  LCL 

CHARACTER* 10  NFILE 

DATA  A2/0.0.  1.88.  1. 02.  .  73,  .  98.  .  48.  .  42. .  37. .  34.  .  31/ 

DATA  Al/O.O,  3.  76.  2.39.  1.88.  1.60.  1.41.  1.28.  1.  17.  1.09.  1.03/ 
DATA  04/0.0.3.27.2.97.2.28.2.11.2.00.  1.92.  1.86.  1.82.  1.78/ 
DATA  D3/0.  O.  0.  O.  O.  O.  O.  O.  0.  0.  O.  O.  .  08,  .  14.  .  18.  .  22/ 

DATA  B4/0.0.  3.27.2.  97.2.27.  2.  09.  1. 97.  1. 88,  1. 81 .  1.76.  1.72/ 
DATA  B3/0.0.  O.  0.  0.0.  0.0.  0.  0.  .03.  .  12.  .  19.  .24. .  28/ 

DATA  D2/0.0.  3.69.  4.36.  4.  70.  4.92.  9.08.  9.  2.  9.  31.  9.  39,  9.47/ 
DATA  D1  /0. 0.  0. 0.  0. 0.  O.  O.  0.  0.  0.  0.  .  2.  .  39.  .  99.  .  69/ 

DATA  SD2/0.  0.  1.  128,  I.  693,  2.  099.  2.  326.  2.  934.  2.  704.  2.  847. 

*  2.  970.  3.  078/ 

DATA  B2/0.0.  1.84.  1.86.  1.81.  1.76.  1.71.  1.67.  1.64.  1.61.  1.98/ 
DATA  B  1/0.0.  0.0.  0.0.  0.0.  0.0.  .03.  .  10.  .  17.  .22.  .26/ 

DATA  C2/0. 0.  .  9642.  .  7236.  .  7979.  .  8407.  .  8686.  .  8882.  .  9027. 

«  .9139.  .9227/ 

WRITE!9.  •>  ’ 


MENU 


c 

c*** 

c*** 

c 


WRITEO.*)  '!>  CONTROL  LIMITS  CALCULATED  FROM  DATA1 
WRITEO.  •>  ' 2)  USE  OLD  CONTROL  LIMITS' 

WRITEO. «)  '3>  USE  AIMED  AT  VALUES' 

READ  <3.0  MENCL 
IF<MENCL.EO. i)  THEN 

IF<MENU2.E0.  3)  THEN 

I F  <  MENU 1 . EO. 1 1  THEN 
UCL-KBB+A2(SS)  *SB 
LCL«XBB-A2(SS)*SB 
CL«XBB 

ELSE  IF (MENU1 • EQ. 2 )  THEN 
Y«(NS-1.  )  /NS 

UCL«XBB*SORT(Y)*Al(SS)*SB 
LCL-XBB-S0RT(Y)*A1 ( SS ) *SB 
CL-X08 
END  IF 

ELSE  ’  r  <  MENU2.  EO.  4.  AND.  MENU  1 .  EO.  1 )  THEN 
:L-D4(SS)«SB 
LCL«D3(SS)*SB 
CL-SB 

ELSE IF < MENU2.  EO.  4.  AND.  MENU  t . EO.  2)  THEN 
UCL«B4(SS)*SB 
LCL-B3(SS)«SB 
CL-SB 
ENDIF 

ELSE IF < MENCL. EQ.  2 )  THEN 
IF  <MENU2. EQ. 3 )  THEN 

OPEN (UN IT-3.  FILE-NFILE. STATUS- 'OLD' > 

READ<3.  •>  UCL.  LCL.  CL 
CLOSE  <UNI T-3. STATUS- ' KEEP ' > 

ELSE IF ( MENU2. EO. 4  >  THEN 

0PEN<UNIT«4. FILE-NFILE. STATUS- 'OLD ' » 

READ<4.  •)  UCL.  LCL.  CL 
CLOSE (UNIT— 4.  STATUS- 'KEEP ' > 

ENDIF 

ELSE IF (MENCL. EO. 3)  THEN 
IF (MENU2. EQ. 3)  THEN 

WRITEO.  »)  'ENTER  AIMED  AT  VALUES  FOR  XBAR  AND 

1  5I0MA ' 

READO.  •)  XB.SIC 
V-REAL(NSI 

UCL-XB+ ( 3/SORT ( Y ) >»SIO 
LCL— XB- ( 3/SQRT  <  Y  » I *SIC 
CL-XB 

ELSEIF(MENU2.E0.4>  THEN 

WRITEO,*)  'ENTER  AIMEO  AT  VALUE  OF  SICHA' 
READO.  O  SIO 
IF (MENU 1 . EO. I >  THEN 
UCL-D2(SS)*SIC 
LCL-D1 (SS).SIG 
CL-SD2(SS)*Sie 
ELSE 

Y— (NS— 1.  ) /NS 

UCL-B2 ( SS  >  *S IQOSQRT ( Y I 

LCL— B 1 <  SS  >  »S 1 0*SORT ( Y ) 

CL— C2(SS I *SIC*SQRT  <  Y) 

ENDIF 

ENDIF 

ENDIF 

RETURN 

END 

THIS  SUBROUTINE  CALCULATES  THE  CONTROL  LIMITS  FOR  A 
P  (AVERASE  LOT  SUE).  NP.  OR  C  CONTROL  CHART 

SUBROUTINE  CLPNPC (MENU I. MNP.  AP. AN. UCL.  LCL.  CL. NFILE) 
REAL  LCL 

CHAR AC TER *10  NFILE 
WRITEO.*)  '  MENU' 

WRITEO.  •>  '  ' 

WRITE  (3.*)  'O  CONTROL  CHARTS  CALCULATED  FROM  DATA ' 
WRITEO.*)  '2>  USE  OLD  CONTROL  LIHTIS' 

WRITEO.*)  '3)  USE  AIMED  AT  VALUES' 

READO.*)  MENCL 
IF(MENCL.EO. 3)  THEN 

IF  <  MENU1 .  EO.  4  >  THEN 

WRITE (3.  *)  'ENTER  AIMED  AT  VALUE  OF  P-PRIME' 
ELSE IF (MENU I . EO.  3)  THEN 

WRITEO*  •>  'ENTER  AIMED  AT  VALUE  OF  C-PRIME' 
ENDIF 

REAOO.  #>  AP 
ENDIF 

IFfMENCL.EO.  1  *  OR. MENCL. EQ. 3>  THEN 
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I F ( MENU 1 . EO. 4 )  THEN 
IF(MNP. EG.  t  )  THEN 

UCL-AP*3*S0RT(  (AP*(  l-AP) >/AN) 

LCL-AP-3»SORT ( <  AP* ! | -AP > > /AN  ) 

CL-AP 

ELSE 

UCL-AN*3*S0RT ( AN* ( I -AP ) ) 

LCL-AN-3*SQRT ( AN* ( 1 -AP  ) ) 

CL- AN 
END  IF 

ELSE IF (MENU!.  EQ.  3>  THEN 
UCL«AP*3*S0RT ( AP ) 

LCL«AP-3*SQRT<AP ) 

CL-AP 
END  IF 

ELSE 

0PEN(UNIT«3.  FJLE-NFILE.  S'!  US- 'OLD') 

READ  ( 3.  * )  UCL.  LCL.  CL 
CLOSE (UN IT-3. STATUS- 'KEEP  , 

END  IF 

RETURN 

END 

THIS  SUBROUTINE.  CALCULATES  THE  CONTROL  LIMITS  FOR  A 
P  (INDIVIDUAL  LOT  SIZES)  OR  U  CHART 

SUBROUTINE  CLPU( MENU I .  N.  NV,  AP.  VUCL.  VLCL.  CL> 

DIMENSION  VUCLOO).  VLCL (30) 

REAL  NVI30) 

WRITEO.  *)  '  MENU' 

WRITE (3.  •)  '  ' 

WRITEO.  *>  'I)  CONTROL  CHARTS  CALCULATED  FROM  DATA' 
WRITEO.*)  '2>  USE  AIMED  AT  VALUE  FOR  CENTRAL  LINE' 

READO.  •)  MENCL 
IF (MENCL. EQ. 2)  THEN 

IF ( MENU 1. EO. 3)  THEN 

WRITEO.  *)  'ENTER  AIMED  AT  VALUE  FOR  U-PRIME' 

ELSE 

WRITE! 9*  *)  'ENTER  AIMED  AT  VALUE  FOR  P-PRIME' 

END  IF 

READO.*)  AP 
ENDIF 

IF (MENU1. EO. 3)  THEN 
DO  10  I-l.N 

VUCL  <  I  > -AP*3*SGRT  (  AP/NW  2  )  ) 

VLCL ( I ) — AP— 3*S0RT ( AP/NV ( I ) ) 

IF ( VLCL ( I ) . LT. 0.0)  VLCL(l)-0.0 
CONTINUE 
CL-AP 

ELSE 

DO  20  I-I.N 

VUCL ( I ) -AP*3*SORT ( AP* ( l -AP ) /NV< I ) ) 

VLCL ( I ) -AP-3*SORT  <  AP* ( l -AP ) /NV ( I ) ) 

IF ( VLCL ( I ) • LT • 0.  O)  VLCL(I)-0.0 
CONTINUE 
CL-AP 
ENDIF 
RETURN 
END 

THIS  SUBROUTINE  PRINTS  X-BAR,  R,  8.  P  (AVERAGE  LOT  SIZES).  NF 
AND  C  CONTROL  CHARTS  TO  THE  SCREEN  OR  TO  AN  OUTPUT  FILE 

SUBROUTINE  CHART! A. UCL.  LCL. CL, N. NFILE.  MENU l .  MENU 2.  MNP.  MCU. •) 
DIMENSION  A (30) .  NBOO) 

REAL  LCL 
INTEOER  86 
CHARACTER* 10  NFILE 
CHARACTER*70  LINE 
CHARACTER* 1  POINT 
WRITEO.*)  '  MENU' 

WRITEO.  •)  '  ' 

WRITE! 9.*)  '1)  CONTROL  CHART  TO  FILE  * 

WRITE! 9. *)  '2)  CONTROL  CHART  TO  SCREEN' 

WRITE (9. •)  '3)  BOTH' 

READ <9,  *)  MENCH 
IF(MENCH.EQ.2>  00  TO  129 

IF (  (MENUS.  EO.  I.  OR.  MENUS.  EO.  2) .  AND.  MENU2.  EO.  3)  THEN 
WRITE! 2.  10) 

FORMAT! 'I '.  39X.  'X-BAR  CHART') 

ELSE IF! MENU I.EO.  I.  AND. MENU2. EO. 4)  THEN 
WRITE! 2.  20) 

FORMAT!  '1 '. 39X.  'RANGE  CHART') 


ELSEIF!MENU1.E0.2. AND.MENU2.EQ.4)  THEN 
URITEO.  30) 

FORMAT!  *1  *.37X.  'S  CHART*) 

ELSE1FCMENU1.EQ.4. ANO.hNP.EQ. I )  THEN 
URITEO. 40) 

FORMAT!  *1  37X.  *P  CHART*) 

ELSE IF (MENU1.EQ. 4. ANO.HNP.EQ. 2)  THEN 
URITEO. 30) 

FORMAT C  *  1  * , 37X.  'NP  CHART*) 

£LSEZF<MENU1.E0.3.  ANO.MCU.EO.  1  >  THEN 
WRITE < 2.  60) 

FORMAT  C  *  I  *  #  37  X .  *C  CHART*) 

ENDIF 

IFILCL.LT.O.  O)  THEN 
WRITE  <  2. 70) 

FORMAT <  *0*.  64X.  *U*/40X.  *C*.24X.  *C'/40X.  'L'.24*. 

I  *L  */) 

URITEO.  731  CL.  UCL 

FORMAT (  //36X.  F9.  4.  16X.  F9.  4/40X.  *;  *.24X.  *:  *) 

ELSE 

mRITEO.  60) 

FORMAT!  'O'.  14X.  *L*.  49X.  *U'/13X.  *C'.24X.  *C*.24X,  *C'/13X. 
I  24X.  *L'.  24X.  'L*/> 

wR  I  TEC  2.  90)  LCL.  CL.UCL 

FORMAT <  1 1 X.  F9.  4.  16X.  F9.  4,  16X.  F9.  4/1 3X.  *:  *.24X.  *:  *.  24X.  ' 
ENDIF 

WRITE (2* 93) 

FORMAT ( 3X.  72!  *.  *)  ) 

DO  113  I-l.N 
L I NE  ■  *  * 

XF<LCL. LT. 0. O)  00  TO  96 
LINEC10: 10)-*! * 

LINE (33: 33)-' ! * 

LIN£(60: 60)-*! ' 

URITEO. 100)  LINE 
FORMAT  <  3X.  *.  A70.  *  *> 

LINE! 10:  10)-'  * 

LINE (60:  60)  - '  * 

J-  <  <  49*A  <  I  )  ♦UCL— SO*LCL  )  /  (  UCL-LCL  )  )  ♦.  3 
J— J*9 

IF!J.LT.  1)  J-l 
IF < J. OT. 70)  J— 70 
XFCJ.LT.  10.  OR.  J.GT.60)  THEN 
LINE! J:  J)- *? ' 

ELSE 

LINE(J:  J)- '• ' 

ENDIF 

WRITE (2. 110)  I*  LINE 

FORMAT (2X.  12.  *.  A70.  '  '> 

CONTINUE 
WRITE (2.  120) 

FORMAT ( 3X.  72(  '.  ') ) 

IF (MENCH. EG. 1 )  CO  TO  260 

IF < ( MENU1 . EQ. 1.0R.MENU1.EQ. 2). AND. MENU2. EO. 3)  THEN 
WRITE! 3.*)  *  X— BAR  CHART' 

ELSEIF(MENU1.EQ. 1. AND.MENU2.EQ. 4)  THEN 
URITEO.*)  *  RANGE  CHART* 

ELSEIF(MENU1.EQ.2.  AND.MENU2. EQ. 4 )  THEN 
URITEO.*)  *  S  CHART* 

ELSEIF (MENU1 . EQ. 4. AND. MNP. EQ. 1 )  THEN 
WRITE <3. •  )  *  P  CHART' 

ELSEIF!MENU1. EQ. 4. AND. MNP. CO. 2)  THEN 
URITEO. *)  *  NP  CHART' 

ELSEIF (MENU1  •  EQ.  3.  AND.  MCU.  EQ.  1 )  THEN 
URITEO.*)  *  C  CHART* 

ENDIF 

DO  210  t-l.N 

NB<  I  )•(  1 16*A(  X  )*UCL-17*LCL > / (UCL-LCL )  )♦.  3 
NB  < I ) -NB ( I > ♦ 1 
IF(NB( I >.LT.  1 )  NB < 1 )  — l 
IF(N8 ( I ).  CT.  19)  NB <  X )  —  1 9 
CONTINUE 
DO  230  1-19.  1.  -I 
DO  213  J-t. 70 

LINE! J: J)*  *  * 

CONTINUE 
IFCI.EQ. 10)  THEN 
DO  220  J-l.  70 

LINE! J:  J)-*- ' 

CONTINUE 

ELSEIF!  I.EQ.2.0R.  I.EQ.  16)  THEN 

IF!  I.EQ.2.  AND.  LCL.LT.O.  0)  CO  TO  227 
DO  223  J-l.  70.  2 


WvV\'-v  *v  \ -\  v.v/'.v.v;.*;.  ** *.**•  .•  v 


WV’ 


v  v'.-j 


L INE!  J  :  J)-'-' 

CONTINUE 
ENOKF 
POINT-'*  ' 

IFU.EO.  19.  OR.  I.  EQ.  I)  POINT*  ' ? 4 
DO  230  J— 1#  N 

IFU.EO.  NB!J)  )  LINE !  J*2:  J*2> -POINT 
CONTINUE 
IFU.EO.  IB)  THEN 

WRITE (3. 240)  UCL.  LINE 
FORMAT C  *  *.F8.4.  *  '.AA1> 

ELSE  IFU.EO.  10)  THEN 

WRITE!  9*  240)  CL.  LINE 
ELSE IF! I. EO. 21  THEN 

IFILCL.LT. 0.0)  GO  TO  290 
WRITE (9. 240)  LCL.  LINE 

ELSE 

WRITE! 9. 249)  LINE 
FORMAT!'  *  *  10X#  A61  > 

END  IF 
CONTINUE 

WRITE!9. *>  'ENTER  1  TO  SAVE  THESE  CONTROL  LIMITS' 
WRITE!9.  *)  'ENTER  2  TO  CONTINUE' 

READ! 9.  *)  NSCL 
IFINSCL. EO. 1 )  THEN 

IF ! ! MENU I. EQ. 1 . OR. MENU2. EO. 2 ) . AND.  MENU2. EQ.  4 >  THEN 
NU— 4 

ELSE 

NU— 3 
ENDIF 

OP£N(UNIT-NU.  FILE-NFILE.  STATUS- 'UNKNOWN ' ) 

CLOSE! UN ZT-NU.  STATUS- 'DELETE ' ) 

OPEN! UNI T-NU# FILE-NFILE.  STATUS- 'NEW ' > 

WRITE!NU.  *)  UCL*  LCL.  CL 
CLOSE! UN I T-NU.  STATUS- 'KEEP ' > 

ENDIF 
RETURN  1 
END 

THIS  SUBROUTINE  PRINTS  P  ! INDIVIDUAL  LOT  SIZES)  OR  U 
CONTROL  CHARTS 

SUBROUTINE  UCHART !A.  VUCL.  VLCL.  CL*  N.  MENU I .  *> 

DIMENSION  A! 30).  VUCL! 30).  VLCL! 30), NA!30>*  NB!30)«  NC!30> 

CHARACTER* 10  NFILE 

CHARACTER *70  LINE 

CHARACTER* 1  POINT 

WRITE!9. •)  '  MENU' 

WRITE! 9. *)  '  ' 

WRITE! 9.*)  'i)  CONTROL  CHART  TO  FILE' 

WRITE! 9. *>  '2>  CONTROL  CHART  TO  SCREEN' 

WRITE! 9.*)  #3>  BOTH' 

READ!  9.*)  MENCH 
XHI-VUCL! 1 ) 

XLO— VLCL! I ) 

DO  9  I-l.N 

IF ! VUCL! I ) • OT . XHI )  XHI-VUCL! I ) 

IF! VLCL! I ). LT. XLO)  XLO-VLCL!I) 

CONTINUE 

IF!MENCH.EQ.2>  CO  TO  90 
IF!MENUi. EQ. 3)  THEN 
WRZTE!2.  10) 

FORMAT!  '1  37X.  'U  CHART') 

ELSE 

WRITE! 2.  20) 

FORMAT!  '1  '.  37X.  'P  CHART') 

ENDIF 

WRITE!2«  30)  CL 

FORMAT!  'O'.  33X.  'CL  -  '.F9.4.//) 

WRITE! 2.  40) 

FORMAT ! 9X» 72!  '.  ' ) ) 

J- ! ! 49*CL*XH I -90*XL0 ) / ! XHI-XLO )  )  ♦.  9 
J-J+9 

DO  70  I-l.N 
LINE-'  ' 

K— ( ! 49* A! I ) ♦XHI- 90* XLO) / ! XHI-XLO ) >♦.  9 
K«K*9 

L«< ! 49*VUCL ! I ) ♦XHI -90* XLO)/ ! XHI-XLO) )♦. 9 
L-L+9 

M-! !49*VLCL! I >*XHI-90*XL0) / ! XHI -XLO) )♦. 9 
M»M*9 

LINE!  J:  U)-'  !  ' 

L INE !L  L )*  '  ‘  ' 


I 


m 


LINE(M.N)-'!  * 

IFIK.LT.  I)  K-l 

IFIH. CT. 70)  K«70 
IFIK.LT. M.OR.K.CT.LI  THEN 

LINEIK:K)-‘T' 

ELSE 

LINEIH.R)-'*' 

END  IF 

URITEO.  90)  I.  LINE 
FQRMATI2X.  12.  '  * .  A70.  ') 

LINEIK: K)- *  ' 

IFIK.EO.  U.  OR.  K.  EO.  L.  OR.  H.  EO.  N>  LINEIH.  K>-'  <  ' 

WHITE <2. 60)  LINE 
FORMAT  OX.  *:  '.  A70.  '  ' ) 

CONTINUE 
URITEO.  SO) 

FORMAT  OX.  72<  '.'ll 
IF (WENCH. EO. 1 )  CO  TO  160 
IF < MENU  1. EO. 3)  THEN 

URITEO.  •>  '  U  CHART' 

ELSE 

URITEO.  •>  '  P  CHART' 

ENDIF 

N0-<  < 16*CL*XHI-17*XL0)/( XHI-XLO) >  +  l. 5 
00  110  I-l.N 

NA( I >-( ( 16*A< I )+XHI — 1 7*XL0 ) / ( XHI— XLO) )♦!. 9 
NB t I >-! < 16*VUCL< I >*XHI-I7*XL0>/ ( XHI-XLO) )♦!. 3 
NCI I)-l  < 1 6* VLCL 1 1 >*XHI-17*XL0)/< XHI-XLO) )♦!. 9 
IF1NA(I).LT.  1 )  NA(I>-1 
IFCNAID.QT.  IV)  NA1D-19 
CONTINUE 
DO  190  I-1V.  l.-l 
DO  119  J-1.70 

LINE!  J:  %!)-  '  ' 

CONTINUE 

IFII. EO.NO)  THEN 
DO  117  J-1.70 

LINE! J: 

CONTINUE 

ENDIF 

DO  120  U-l.N 

IF(  I.EO.  NB(  J).  OR.  I.  EO.  NCI  J)  )  THEN 
LINEIJ02:  J*2>-'-' 

LINE! J»2+l : J»2+l )» ' 

ENDIF 

POINT-'*' 

IF1NA1 J)  ■  LT.  NCI J).  OR. NA( J). CT. NB ( J) )  POINT-'?' 
IF! I. EO. NAl J) )  LINE! J*2+l : J*2+l )— POINT 
CONTINUE 
IF1I.E0.ND)  THEN 

URITEO.  130)  CL. LINE 
FORMAT!'  '.FV.4.  '.AAD 

ELSE 

URITEO.  140)  LINE 
FORMAT!'  ’.  10X.  A61) 

ENDIF 

CONTINUE 

URITEO.  * )  'ENTER  1  TO  CONTINUE' 

READI9.  •>  NK 
RETURN  1 
END 


Setts' 
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